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LOOKING FOR SINGLE TOP QUARKS IN D0 DATA 
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The search for the electroweak production of top quarks is entering an exciting 
time: with the ever increasing luminosity furnished by the Tevatron pp CoUider at 
Fermilab and sophisticated analysis techniques we are now able to approach the 
expected standard model production cross section for this new mode of production. 
Using 370 pb~^ of D0 data and likelihoods as discriminants to extract the signal 
from the large backgrounds, we set upper limits on the standard model production 
cross section in the s-channel of 5.0 pb and in the t-channel of 4.4 pb, at 95% 
confidence level. 

1 Introduction 

Top quarks have been observed so far produced in tt pairs in collisions at 
Fermilab's proton-antiproton Tevatron Collider. Pair production occurs pre- 
dominantly when two incoming quarks exchange a gluon {qq ^ g ^ tt): that 
is, they are mediated by the strong force. But another mode of production, 
which has not been detected yet, is predicted by the standard model (SM), 
when the incoming quarks exchange a Vl^-boson and as a result a single top 
quark is produced along with a 6-quark or other light quarks. This is an elec- 
troweak top quark production, also referred as single top quark production. 
There are two main modes of W boson exchange with top quark production 
observable at the Tevatron: s-channel and t-channel, as depicted in Fig. [TJ 
At NNLO, the predicted cross sections are 0.88 and 1.98 pb respectively. 




Figure 1. The two dominant production modes of single top quarks at the Tevatron: s- 
channel (left) and t-channel (right) exchange of a boson. 



The observation and study of this new production mechanism is impor- 
tant for the foUowing reasons: a measurement of the cross section can yield 
the first direct measurement of the CKM matrix element |Vt6| without any as- 
sumptions on the number of families in the SM; new particles can enhance the 
production cross section in the s-channel, and anomalous couplings can alter 
the production through the t-channel, thus this channel offers a very sensitive 
window into new physics; this channel can be used to test the V-A nature of 
the SM; and not least important, this channel requires a good knowledge of 
the W-|-jets background and is itself a background to Higgs boson searches. 
For more details of the analyses described here, see Ref. 

2 Signal and Background Topologies 

The signal topology we are interested in arises from the top quark decaying to 
Wb and the W decaying leptonically to electrons or muons. We are looking 
for one isolated high lepton, missing transverse energy, and at least one 
6-quark jet (two in the case of the s-channel, and usually only one in the 
t-channel, since the b from gluon splitting tends to get lost in the beam-pipe). 
The main background is W-|-jets production: these events tend to be a little 
bit more pencil-like and tend to have smaller total reconstructed energy. But 
the cross section for W-|-jets is at least two orders of magnitude bigger than 
that of single top. The second largest background is ti production in the -^-l-jets 
or dilepton decay modes. These events tend to have more total reconstructed 
energy (having two heavy top quarks in them) and be more spherical than 
single top. The ti cross section is around twice that of single top. 

The iy-|-jets background is modeled using Alpgen^^, but it is normalized 
to the data. The tt background is also modeled using Alpgen and normalized 
to the NNLO cross section of 6.7 pb. The QCD background is obtained 
from data and the fraction present in the signal data is obtained from the 
same normalization of W-|-jets to signal data. For single top, we employ 
CompHEP ^, which effectively describes all angular correlations between the 
final state particles at NLO. 

3 Analyzed Data and Event Selection 

Data collected with the D0 detector ^ between August 2002 and October 
2004 corresponding to nearly 370 pb~^ have been analyzed in the search for 
s- and i-channel single top quark signals. 

We apply loose selection cuts to the data and the Monte Carlo to extract 



signal-like physics: we require to have 2, 3 or 4 jets in the event, the leading one 
with a pt > 25 GeV and < 2.5, the rest with pr > 15 GeV and |ry| < 3.4; 
we additionally require one isolated electron or muon with > 15 GeV 
and |77(e)| < 1.1 or |77(/i|) < 2.0; and the transverse missing energy to be 
> 15 GeV. 

Finally, we require the presence of at least one jet identified as coming 
from a b quark. We use the Jet Lifetime Probability (JLIP) algorithm based 
on the fact that jets originated from a b quark exhibit tracks with large impact 
parameters in the transverse plane and thus allow to separate jets originated 
close to the primary vertex or not. For each jet we have a continuous distri- 
bution and decide to "tag" the jet as a tight b jet if Pjlip < 0.3% or loose b 
jet if Pjlip < 1%- The samples are separated into events with one tight and 
no loose 6 jets, and events with one tight b jet and one loose b jet. Allowing 
the second b tagged jet to be looser than the first, increases our acceptance 
and somewhat reduces the uncertainty coming from b tagging in these small 
samples. 



4 Discriminant Variables and Likelihood discriminants 

At this point in the analysis, cross-check samples enriched in VF-|-jets and tt 
events only are utilized to convince ourselves that the background modeling is 
correct not only in the signal region but also in other regions. Figure [2] shows 
some of the variables in the signal region used later in the analysis. Good 
agreement is found in all of them between data and background modeling. 

Single top is a difficult signal to extract from the large backgrounds, as 
can be seen from Fig. [21 Since no single distribution allows to separate the 
signals from the backgrounds, multivariate analysis techniques are best suited 
to efficiently examine a large portion of the phase space trying to identify and 
isolate regions were the signal sticks out. For this analysis we have employed 
likelihood discriminants with mostly uncorrelated variables as input. Here is 
where we benefit from the loose selection cuts: the likelihood discriminants do 
a better job of extracting the signal than any optimization of simple stringent 
cuts on basic object distributions. 

The likelihood discriminant is defined as: C(x) = ~ ^-T^'g'"'''^'' t^y, 

' signalV^j"r' background j 

where 7-*signai(a;) and 7-'background(a;) are the probability density functions of 
signal and background events as derived from the product of input template 
distributions. We have performed this analysis in separate channels for elec- 
tron/muon, single/double tagged, s-channel/f-channel signals and W+iets/tt 
backgrounds giving a total of 16 different likelihoods. Given the different 
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Figure 2. Distributions for data and background with the electron and muon channels 
combined of: leading jet pxi the ^x-l-lepton px, the charge of the lepton times the eta of 
the untagged jet, and the invariant mass of the reconstructed W and the leading tagged 
jet. 



topologies of the main backgrounds and signals, we thus benefit from the dif- 
ferent separation power that each likelihood can accomplish in each channel 
and then can easily combine them for the extraction of the limit. The output 
of some of the likelihoods can be seen in Fig. [31 

5 Results and conclusions 

We use two-dimensional distributions of the output likelihoods of tt vs. 
W+iets for each signal, and use the shape information on each bin of that 
plane, taking into account the correlations in the errors in each bin, in a 
Bayesian analysis with a Poisson distribution for the observed counts in data 
and flat priors for the signal cross sections. The priors for the signal ac- 
ceptance and the backgrounds are multivariate gaussians. The observed (ex- 
pected) 95% confidence level upper limits are 5.0 (3.3) pb for the s-channel 
and 4.4 (4.3) pb for the t-channel. 

With more than 1 fb~^ of delivered luminosity and better object recon- 
struction, calorimeter calibration, much reduced jet energy scale errors (which 
improve our invariant mass calculation), better ^-tagging techniques, and more 
powerful separation techniques like matrix element and boosted decision trees, 
we can expect the single top search at D0 to enter a very exciting time in the 
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Figure 3. Likelihood outputs for electron and muon channels combined in the single tagged 
tqb/tt (upper left) and tqb/W+]ets (upper right) filters; and for double tagged tb/tt (lower 
left) and tb/W +}cts (lower right) filters. 



near future. We are at the same time looking for physics beyond the standard 
model in this final state topology coming from a new massive W boson and 
flavor changing neutral currents-^. 
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